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Noiseless non-reciprocity in a parametric 
active device 


Archana Kamal 1 , John Clarke 2 and M. H. Devoret 1 * * 

Non-reciprocal devices such as circulators and isolators belong 
to an important class of microwave components employed in 
applications including the measurement of mesoscopic circuits 
at cryogenic temperatures 1-5 . The measurement protocols 
usually involve an amplification chain that relies on circulators 
to separate input and output channels and to suppress 
backaction from different stages on the sample under test. 

In these devices the usual reciprocal symmetry of circuits 
is broken by the phenomenon of Faraday rotation based on 
magnetic materials and fields 6 . However, magnets are averse 
to on-chip integration, and magnetic fields are deleterious to 
delicate superconducting devices 7 ' 8 . Here we present a new 
proposal that combines two stages of parametric modulation 
to emulate the action of a circulator. It is devoid of magnetic 
components and suitable for on-chip integration. As the 
design is free of any dissipative elements and based on 
reversible operation, the device operates noiselessly, giving 
it an important advantage over other non-reciprocal active 
devices for quantum information processing applications. 

Reciprocity is one of the fundamental symmetries frequently 
encountered in electrical circuits. It is equivalent to the more 
familiar notion of the principle of reversibility in optics, which 
states that any experiment is symmetric under an exchange of 
source and image 9 . Reciprocity can, however, be violated, for 
example, by the magneto-optic Faraday effect 6 —the rotation of the 
polarization vector of light resulting from different propagation 
velocities of left- and right-circularly polarized waves in the 
presence of an applied magnetic field B parallel to the direction of 
propagation (Fig. 1). The non-reciprocal phenomenon of Faraday 
rotation should be contrasted with the superficially similar, though 
reciprocal, effect of optical activity where the polarization vector of 
light is rotated on passage through a non-centrosymmetric (chiral) 
medium. This change in the sense of rotation of polarization for 
counterpropagating waves in a Faraday medium (as seen by the 
observer receiving the light) has led some physicists to refer to this 
effect as a form of time-reversal symmetry breaking, a use of words 
that we prefer to avoid here 10 . 

The phenomenon of non-reciprocity has propelled numerous 
theoretical investigations (ref. 11 and references therein); further¬ 
more, it offers immediate practical applications. Recent progress 
in solid-state superconducting qubits, which provide some of the 
most promising architectures for scalable quantum computers 12 ’ 13 , 
has generated a strong incentive to integrate the components 
required for qubit operations and readout on-chip for incorpora¬ 
tion in future quantum mechanical processors. A large variety of 
qubit readout protocols involve microwave-reflection-based mea¬ 
surements and rely on non-reciprocal devices such as circulators 
(or isolators) for separation of input and output channels 3,5,14 . 
These devices also play a strategic role in measurements based 




Figure 1 | Faraday rotation and circulator action, a, Faraday rotation for a 
wave travelling from left to right in a Faraday-active medium, followed by a 
reflection back into the medium leading to a reversal of the direction of 
propagation. The rotation of the light polarization is fixed to a rotation-like 
property of the medium (shown by the arrows), set by an external magnetic 
field oriented along the propagation axis. The sense of light rotation as seen 
with respect to the direction of propagation remains the same, leading to 
the doubling of the rotation angle on reversing the ray through the medium, 
b, Rotation of the polarization vector of light on passage through an 
optically active medium, on the other hand, cancels out on reversing the 
direction of propagation. This occurs because optical rotation depends on 
the chirality of the medium (represented as a helix), which also reverses 
with the direction of propagation, c, Representation and schematic design 
of a conventional four-port circulator. The device consists of two 90° 
hybrids (equivalent to optical beam splitters) separated by a non-reciprocal 
phase shifter based on Faraday rotation. The filled black arrows indicate an 
amplitude split with no phase change and the open arrows indicate an 
amplitude split with a 90° phase change. The non-reciprocal phase shift is 
effective only for the propagation direction indicated by the arrow on the 
phase-shifter rectangle. 
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Figure 2 | Description of an active reversible (information-conserving) 
in-phase/quadrature modulator carrying out frequency UDC. a, Circuit 
schematic of the UDC containing only dispersive components. The two 
low-frequency series LC resonators (with L/\ = Lb and Ca = Cb) are fed by 
two input semi-infinite transmission lines, A and B, and parametrically 
coupled to a third high-frequency series LC resonator leading to an output 
line C. The parametric coupling is achieved by varying the mutual 
inductances Mi and M 2 between the left and right resonators at the carrier 
frequency co c , which, for optimal frequency conversion, is set at the band 
centre of the right resonator. When operated from left to right, the circuit 
carries out the modulation of low-frequency signals of frequency u>q 
travelling on ports A and B to generate sidebands at w c ±a >0 travelling on 
the high-frequency line C. It carries out the inverse operation of 
demodulation when operated in reverse from right to left, b, Spectral 
density/response landscape for different spatial channels (or ports) of the 
UDC circuit in a as a function of frequency. The dotted lines represent the 
couplings between different ports. The solid and the dashed arrows 
represent different frequencies and respective conjugates. The resonance 
line shapes of the two spatially distinct ports A and B are centred at 
a > a = o)b = 1 AAa.bCa.b- ^ ere we s how the case when the incoming signal 
at a>o is resonant with the centre frequency ( 0)0 =o>a.b)- The two 
sidebands generated by the UDC on channel C are detuned from the pump 
at u) c by equal amounts. 


The main building block of our design is a reversible 
in-phase/quadrature modulator capable of carrying out noise¬ 
less frequency up- and down-conversion. A convenient analyt¬ 
ical model capturing the fundamental properties of the device 
is shown in Fig. 2a. The device comprises two low-frequency 
inductor-capacitor (LC) resonators (addressed by two semi-infinite 
transmission lines A and B) coupled to a high-frequency resonator 
(addressed by the transmission line C) through time-varying cou¬ 
plings Mi,M 2 that emulate the role of the pump drive in active 
nonlinear devices and transfer energy from the tone at a) c to the 
signal modes propagating on the transmission lines. It operates 
in a manner analogous to the in-phase/quadrature modulation 
schemes routinely used in radiofrequency communication systems 
and microwave pulse engineering (hence the name) and converts 
two orthogonal spatial modes travelling on two distinct spatial ports 
(A, B) at the same frequency (here 000 ) into two orthogonal temporal 
modes travelling on the same spatial line (C) at different frequencies 
(o) + , <»_). However, in view of the reversible and noiseless frequency 
conversion carried out by this device (Fig. 2a), we will henceforth 
refer to it as the up/down-converter (UDC). In practice, such a de¬ 
vice can be implemented on-chip using two ring modulators based 
on Josephson junctions, along the lines of the recently demonstrated 
experiment with a Josephson parametric converter 17 ' 18 . 

The complete design for the active circulator (Fig. 3a) consists of 
a UDC functioning as a frequency up-converter, a phase shifter and 
a second UDC functioning as a frequency down-converter. 

A concise representation of the dynamics at each of the three 
stages in the cascade is provided by the scattering matrix S 
that relates the outgoing wave amplitudes to the incoming wave 
amplitudes as seen from the ports of a network. We begin by 
deriving the scattering matrix of the UDC stage and obtain (see the 
equations (5) and (6) of Methods section) 
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Here a and b denote the (reduced) amplitudes or the annihilation 
operators for the waves travelling on left and right transmission 
lines of the UDC respectively. This notation will be consistently 
used to denote the left- (a) and right- ( b ) propagating modes at 
each stage of the device hereafter (see Fig. 3a). These satisfy bosonic 
commutation relations of the form 19 


on low-noise microwave parametric amplifiers, which, with the 
exception of designs based on the current-biased d.c. supercon¬ 
ducting quantum interference device (SQUID), are also operated 
in a reflection mode with both the input and output signals 
collected on the same spatial channel 15 ' 16 . However, circulators 
(and isolators) routinely use bulk components made from ferrites 
to achieve non-reciprocal phase shifts (Fig. lc) through Faraday 
rotation, making them unamenable to chip fabrication. Moreover, 
to bias the magnetic field in the ferrite, most of these devices use a 
permanent magnet, which may channel flux into the superconduct¬ 
ing device under test. 

Here, we present the full analysis of a model for a four-port 
circulator based on parametric active devices with no magnetic 
components. In active devices the energy source—provided by the 
pump—acts as the external ‘bias’ field and sets the reference phase 
for the system, in analogy with the role played by the magnetic 
field in a Faraday medium. We exploit this effect in a cascade of 
active devices with pump phases at each stage shifted appropriately 
to obtain non-reciprocal transmission. 


[ar /ont , fl/ ln/out ] = S(cOi - a>j) (2) 

In writing equations (1) and (2), we have set a 0 = a[co 0 ],a + = 
= fl[a>_] (see Fig. 2b). Similarly, the reflection coefficients 
at various ports are denoted by r 0 , r + and r_. The cross reflection 
between the low-frequency signal ports is denoted by q 0 . The 
transmission coefficients are written as t (transmission without 
conjugation) and s (transmission with conjugation) with subscripts 
(u, d) indicating the up-conversion and down-conversion, respec¬ 
tively. These scattering coefficients are a function of frequency in 
general (see Supplementary Information for a detailed derivation). 
It is useful to note that the phase of the carrier, denoted by (j> , affects 
only the transmitted amplitudes and rotates the two sidebands in 
opposite directions as can be seen from the corresponding scattering 
coefficients 5 and t in equation (1). The invariance of reflection 
amplitudes to the phase of the coupling will be important in 
understanding total reflections of the cascade, as we describe later. 

Furthermore, we note that the matrix obtained in equation (1) 
is non-unitary, that is S f S / 1, which implies non-conservation 
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Figure 3 | Description of the active circulator, a, Circuit schematic of the active circulator design. The first UDC stage acts as a frequency up-converter 
(UC; also indicated by a gradation in the colour of the relevant square) with a parametric coupling modulated at the carrier frequency 
a>c — ru-|_ — coq = tii— +a>o and a phase 4>l = 0. This is followed by a phase shifter (PS) that phase shifts both the sidebands by ji/2, in opposite directions. 
They are then demodulated by the final UDC stage acting as a frequency down-converter (DC), with the pump phase <p r = it/4, b, Forward (green) and 
backward (red) propagation diagrams calculated using the transfer matrix method for a with an appropriate choice of detuning (<5-t = 1 /VT> and coupling 
strengths («l = “R = 2 -3 / 4 ) for maximum isolation, c, Representation of the device operation using modulation ellipses at each stage in the cascade. The 
top panel shows the forward propagation for the case when two distinct signals enter ports 1 and 2 respectively (a o =1,oo = 0.5e 1 ''; tj = n/24) and the 
bottom panel shows the backward propagation dynamics when the same signals enter ports 3 and 4 respectively (bo = 1,bo = 0.5e"U. The relative phase 
and amplitudes are chosen to represent the most general case of two input signals that differ in both amplitude and phase. The relative phase difference 
between the two signals is encoded as the tilt of the modulation ellipse in the in-phase/quadrature plane and their average phase is represented as the 
colour along the perimeter of the ellipse, with yellow indicating zero phase (also see Supplementary Information for more details). 


of photon number as is natural for an active device. The matrix 
recovers its unitary form as we turn off the couplings Mi,M 2 
responsible for energy transfer between the pump and the signal 
modes. The full 8x8 matrix (Supplementary Equation SI7), de¬ 
scribing the device operation for all modes and their respective con¬ 
jugates, fulfils the fundamental requirement of symplecticity 17,20 . 

We can similarly describe the action of the frequency- 
independent phase-shifting stage using a scattering matrix of 
the form 


/ a + out \ 
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The scattering matrix of the whole device is then obtained using the 
transfer matrix formulation, as described in the Methods section 
(also see Supplementary Information). 

In our analysis we consider the operation at resonance, that is, 
when the input signal frequency coincides with the band centre of 
the input resonators. Setting the phase of the pump at the first UDC 
stage 4 >l = 0 for calculational simplicity, we observe a transmission 
resonance for 8 = ±Jt/2 (phase rotation by the phase shifting stage), 
(j> R = j[/4 (phase of the pump at the second UDC stage), = 1/a/2 
(detuning of the sidebands from the pump in units of linewidth, 
that is, half-width at full-maximum of the resonance line shape) of 
the high-frequency resonator) and a R = a R = M 0 /^/La,bTc = 2~ 3/4 


(strength of the parametric coupling). For this choice of parameters, 
we obtain the scattering matrix of the complete device as 

/ 0 0 0 i \ 

c = 0 0 -i 0 (A) 

itotal i 0 0 0 1 J 

\0 i 0 0/ 

This is the matrix of a perfect four-port circulator. The analogy 
between a conventional circulator and the active circulator design 
proposed in this Letter is made apparent from the respective 
wave propagation diagrams in Figs lc and 3b (see Supplementary 
Information for details on the calculation of coefficients on different 
arms in Fig. 3b). Nonetheless there are important differences 
between the two designs despite the identity of the final S matrix. 
The coefficients on the forward- (green) and backward- (red) 
propagating arms of the active circulator design (Fig. 3b) involve 
deamplification followed by amplification, unlike the passive 
splitters (90° or 180° hybrids) employed in Faraday rotation 
schemes. This can be observed by squaring the amplitudes on each 
of the two arms originating from (or terminating into) a port and 
calculating the net power output, for each isolated UDC stage. It is 
straightforward to observe that, unlike the case of Fig. lc, they do 
not add up to unity. The overall transmission, however, is still unity 
owing to an exact cancellation of the larger and smaller than unity 
gains. The wave propagation diagrams in Fig. 3b reveal another 
important difference of this design from that of a conventional 
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Figure 4 | Variation of the difference between forward and backward transmission coefficients (IS 31 I — |si 3 1 ). a-c, Asymmetry in transmission calculated 
as a function of strength of the coupling a and phase rotation 9 carried out by the second phase-shifting stage (a), detuning <5-t of the sidebands from the 
pump and phase rotation 9 (b), and detuning <5± and the phase of the pump at second UDC stage 0 r (c). The points of maxima correspond to the ideal 
values reported in the text. The plot in b also shows the periodicity of the response of the device as a function of 9. In c, the variation with respect to the 
pump phase shows the reversal of transmission characteristics with 0r ->• —0r. As in b, the response is periodic in 0r with a period equal to jt. It can be 
seen that the design continues to work for moderate deviations from the preferred phase angle 9 = n/2, coupling a\_ = a% = 2 -3 / 4 , detuning = 1/V2 and 
pump phase 0r = jt/4 (values indicated with dashed arrows along the axis). 


circulator. The non-reciprocal action of the active circulator is not 
based on any non-reciprocal phase shifters; instead it relies on the 
active stages used for frequency up- and down-conversion. The 
phase-matching condition in the forward direction is met by tuning 
the phase of the coupling at the input and output UDC stages. In the 
reverse direction, the phase mismatch leads to unity transmission in 
the spatially orthogonal port instead, leading to complete isolation 
between the incident signal port and its corresponding output port. 
Figure 3c shows a convenient method to visualize this circulator 
action geometrically by mapping the device dynamics at different 
stages using a modulation ellipse representation (see the Methods 
section and Supplementary Fig. S2 for details). 

Furthermore, as seen from equation (1), the reflection coeffi¬ 
cients at the UDC stages are non-zero for all modes. Flowever, 
for the whole cascade, the total reflection is identically zero at 
every port (s, t = 0 for all i in equation (4)). This remarkable 
cancellation of total reflections for the cascade can be understood 
in analogy with a Fabry-Perot resonance where a cavity flanked 
by two identical reflecting mirrors exhibits unity transmission at 
resonance. The total phase shift between the active ‘mirrors’ in our 
device: (jc/2) 0+ — (— Ji/2) a t = Jt, is akin to the resonance condition 
when a half-wavelength of the incident radiation equals the length 
of the Fabry-Perot cavity. Also the reflections at the two UDC stages 
are identical (as the reflection coefficients are independent of the 
phase angle <0, equation (1)), fulfilling the second condition for the 
transmission resonance and net cancellation of reflections. A point 
of distinction between the two pictures is that the symmetry of the 
active circulator design is described by a subgroup of the SU( 4) 
group (the group formed by 4 x 4 complex matrices of unit deter¬ 
minant, equation (4)), and not by the SU(2) group that describes 
passive lossless two-port devices such as the Fabry-Perot resonator. 

We note a recent theoretical paper 21 that also showed the 
existence of noiseless non-reciprocity in a nonlinear circuit in¬ 
volving microwave resonators coupled through Josephson junc¬ 
tions. Although this proposal is based on a passive rather than 
active Josephson circuit, unlike our proposal, it involves very small 
Josephson junctions susceptible to offset charges. Also, owing to 
its low characteristic energy, it handles qubit readout signals with 
a lower dynamic range than our proposed device. Furthermore, 
the dependence of the action of our active circulator design on 
different parameters (Fig. 4) shows that the isolation achieved is 
robust to reasonable deviations of parameters from their ideal 
values; thus, the active circulator design holds promise for use in 
practical circuits. Another interesting feature of our device is the 


reversal of its transmission characteristics with the phase of the 
pumps (0 >l.r i-> —<t> l.r) (Fig. 4c). In the classic circulators based on 
passive Faraday rotation, this can be accomplished by changing the 
polarity of the magnetic bias field. Thus, the relative shift of the 
pump phase in an active device indeed plays a role equivalent to 
the magnetic field in a Faraday medium. 

Our scheme for non-reciprocal wave propagation relies on 
the directionality imposed by shifting the phases of the pump 
drives in a cascade of parametrically active stages. The proposed 
design operates noiselessly as the complete device has unity gain 
and consists of purely dispersive components with no dissipation, 
making it attractive for quantum information applications using 
superconducting circuits 22 . As the proposed protocol is generic in 
nature, it should be easily adaptable to optical frequencies where it 
can complement the recently proposed designs of non-reciprocal 
light propagation based on dynamical modulation of the refractive 
index of photonic structures 23 and the use of a surface waveguide 
on photonic crystals 24 . In addition to providing a design for 
practical on-chip non-reciprocal devices, the treatment described 
in this Letter may also give theoretical insights into the inherently 
directional dynamics of devices such as the microwave d.c. SQUID 
(ref. 25), when additional active stages are included in the chain. 
This would be useful in tackling unanswered questions pertaining 
to the quantum noise of d.c. SQUID amplifiers 26 . 

Methods 

The scattering matrix for the UDC is calculated from the network impedance 
matrix Z, as seen from its ports, using the identity 

S = (Z + Z„)- 1 x(Z-Z„) (5) 

where 

Zo = diag(ZA,ZB,Zc,Zc) (6) 

Here, the symbol x denotes matrix multiplication and Za = Zb , Zc are the 
characteristic impedances of the semi-infinite transmission lines serving as low- 
and high-frequency ports respectively. 

For each stage of the cascade, we now transform from the scattering matrix, 
which provides a port-neutral relation between various input and output ports, to 
a port-specific transfer matrix: 



The transfer matrix representation 6 is especially useful for analysing the 
proposed design, because it is straightforward to calculate the total transfer 
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matrix of the device by multiplying the respective transfer matrices of different 
stages in the cascade, 

Total = T DCr X Tps X TjC l (8) 

Here the subscripts L, R index the left-hand upconversion (UC) and right-hand 
downconversion (DC) stage and the subscript PS refers to the intermediate phase 
shifter (Fig. 3a). The scattering matrix of the whole device is obtained from 
Total using the inverse of the transformation in equation (7). For the purpose 
of evaluating Total from equation (8), Tjc:, T>c and T PS are calculated from 
equations (1) and (3) respectively, while noting that 

T>c = F _I x T;c _1 x F 

= PxT oc - 1 xF, (F"‘ = F) (9) 

In equation (9), F = a*® I 2 ( cr x is the two-dimensional Pauli spin matrix and 
I 2 is the two-dimensional unity matrix). This matrix F is required to flip the 
indices, thus maintaining consistency in labelling the ‘in’ and ‘out’ amplitudes 
along a given direction of propagation (see Supplementary Information for 
additional details). 

The modulation ellipses in Fig. 3c are calculated using the transfer matrices, 
with coefficients indicated in Fig. 3b. This representation is inspired by the 
polarization ellipse used to represent the state of polarization of an electromagnetic 
wave (linear, circular or elliptical), which involves recording the trajectory traced 
out by the tip of the polarization vector of light (defined by the instantaneous 
direction of the electric field vector E) in a plane perpendicular to the direction of 
propagation. Equivalently, a two-dimensional representation of the components E x 
and E Y of the electric field in the complex plane can be used to obtain a geometric 
description of the polarization of the light wave. In the case of a modulation ellipse 
representation of Fig. 3 c, we extend this idea to map two distinct orthogonal 
modes (x,y) at each stage of the proposed device (spatial: (x = af ,y = af ) or 
temporal: (x = a+,y = a_)) as an ellipse in the plane defined by the coordinates 
I = Re[x +y], Q= lm[x — y*]. This exercise shows that the final ellipses obtained 
at the output in the cases of forward and backward propagation through the device 
are rotated by 90° with respect to each other, indicating that in the case of reverse 
propagation the orthogonal spatial port, relative to the forward propagation, 
receives the transmitted energy, leading to a circulator action. 
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